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Abstract. We investigate the "low-i^" method (developed by the CCFR/NUTEV collaborations) to deter- 
mine the neutrino flux in a wide band neutrino beam at very low energies, a region of interest to neutrino 
oscillations experiments. Events with low hadronic final state energy v < Vcut (of 1, 2 and 5 GeV) were 
used by the MINOS collaboration to determine the neutrino flux in their measurements of neutrino (//p) 
and antineutrino (f/f^) total cross sections. The lowest i/fj, energy for which the method was used in MINOS 
is 3.5 GeV, and the lowest i?^ energy is 6 GeV. At these energies, the cross sections are dominated by 
inelastic processes. We investigate the application of the method to determine the neutrino flux for i^^, 
energies as low as 0.7 GeV where the cross sections are dominated by quasielastic scattering and Zi(1232) 
resonance production. We flnd that the method can be extended to low energies by using i/cut values of 
0.25 and 0.50 GeV, which is feasible in fully active neutrino detectors such as MINERvA. 
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1 Introduction 

A detailed understanding of neutrino (i/^) and antineu- 
trino {D^) interaction cross sections for various final states 
is required for the next generation neutrino oscillations ex- 
periments. The relevant neutrino energy region of interest 
for the large neutrino detectors such as T2KP, MIN0S[1 
[3], and NOVA 3] is 0.5 < < 3 GeV. 

The MINERvA [S] experiment at the NUMI wide band 
beam at Fermilab uses a fine grain fully active scintillator 
target-detector to investigate neutrino (i^^) and antineu- 
trino (Pf^) cross sections for energies above 0.5 GeV. These 
measurements require a reliable determination of the flux 
as a function of z/^, energy. 

Previous neutrino experiments in wide band beams 
used five methods for the determination of flux as a func- 
tion of energy. 

1. Modeling the distribution of pions and kaons produced 
by incident proton beam in the target. Then, track- 
ing the pions and kaons though the Horn focussing 
magnetic fields, and modeling the decays of pions and 
kaons in the decay pipe. 

2. Measuring the muon fiux that exits the decay pipe and 
relating it to the neutrino fiux. 

3. Monitoring Inverse muon decay events ( j/^ + e — > /i^ + 
Vf,) in the detector. 

4. Monitoring neutrino-electron scattering events (i^^ + 
e ^ Vn + e) in the detector. 



5. The "low-z^" method for the determination of the en- 
ergy dependence of the relative neutrino and antineu- 
trino fiux. 

Here, "low-z^" refers to events with low energy 

transfer to the target nucleon in the scattering processes 
i^f, + N ^ + X and + N ^ fi+ + X. This energy 
transfer manifests itself as the energy {ly = E^ad) of the 
final state hadrons (X) in the laboratory frame. 

There are inherent difficulties in the each of those tech- 
niques: 

1. In method 1, the differential cross sections for the pro- 
duction of pions and kaons by protons incident on a 
thick nuclear target must be known very well. In addi- 
tion, the magnetic field of the horn focussing magnets 
must be modeled reliably. 

2. In method 2, the response of the muon detectors at the 
end of the decay pipe must be very well understood (for 
absolute calibration of the neutrino flux) . The response 
of the muon detectors is sensitive to S rays. In addition, 
since the energy of the muons is not measured, it is 
difficult to determine the energy dependence of the 
neutrino flux. 

3. In method 3, the threshold for the reaction + e — > 
/i^ + i^e is about 12 GeV. Therefore, this method can 
only be used at higher energies. Unfortunately, this 
method cannot be used for the determination of the 
flux for antineutrinos. Inverse muon decay was used 
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by NOMAD to constrain their neutrino flux at high 
energies. 

4. In method 4, only the sum of the fluxes for neutrinos 
and antineutrinos can be measured. This is because 
calorimetric detectors such as MINERvA cannot deter- 
mine the charge of final state electron in J^^+e — )• ly^+e 
events. 

Both methods 4 and 5 are statistically limited. In ad- 
dition, in both methods, the total final state energy in the 
events is not fully reconstructed since there is a neutrino 
in the final state. This places a limitation on the deter- 
mination of the energy dependence of the neutrino and 
antineutrino fluxes. Despite of these limitations, these two 
additional methods are valuable as important consistency 
checks. 

Consequently, having another independent technique 
such as the "low-z/" method is extremely valuable. 

The "low-t"" method was initially developed by the 
CCFR/NUTEV [6l[7] collaboration. At high energy, the 
method was used to determine the relative neutrino flux 
as a function of neutrino energy (E^). The method relies 
on the observation that the charged current differential 
cross section, ^^^r in the limit ly 0, only depends on 
the structure function J-2 , and therefore is independent of 
energy. 

The "low-z^" method was used by the CCFR/NuTeV 
collaborations to measure the energy dependence of a„/E 
and C7i}/E for charged current interactions for energies 
higher than 30 GeV for an iron target. The absolute level 
of the charged cross sections is normalized to previous 
measurements of Uu/ E in a high energy narrow band neu- 
trino beam. 

Most recently, the method was extended to lower ener- 
gies by the MINOS [3] collaboration. The lowest neutrino 
energy for which this method was used in MINOS is 3.5 
GeV for neutrinos and 6 GeV for antineutrinos. 

The absolute normalization used by MINOS is to the 
world average value of charged current Uu/E measure- 
ments for an isoscalar target for neutrino energies between 
30 to 50 GeV. The average value used by MINOS is 



30-50 



= 0.675 X 10-^*cm7GeV 



per nucleon. The antineutrino sample is not independently 
normalized but is related to the neutrinos by using the 
same normalization factor. 

In this communication we investigate the application 
of the technique to much lower neutrino energies {E^ > 
0.5 GeV). Neutrino interactions in this energy range are 
currently being studied at MINERvA. 



2 The "low-z/" method at high energies 

If we neglect terms which are proportional to the niuon 
mass, the differential cross section 



dxdy 



for charged cur- 
rent scattering of (j?^) with an incident energy E^,, 
muon flnal energy E^ and scattering angle 9 can be writ- 
ten in terms of the structure functions Fi = MWi {x, Q^), 
T2 - vm{x, Q2) and = vWz{x, Q^): 



dxdy 



TT 



Mx, 
2E/ 



2 1 I ( 2Mx \2 

2 V l + R 



y 



x:f^), (1) 



where Gp in the Fermi weak coupling constant, M is the 
proton mass, y = v/E^, v = Ehad ^ Ey ~ E^^, Q'^ ^ 
AEyE^sm'^{9/2) is the square of four momentum transfer, 
and X = / {2Mv) is the Bjorken scaling variable. The 
plus sign in front of the xF^, term is for neutrinos and the 
minus is for antineutrinos. 

Here, Q^) is defined as the ratio of the longitudi- 
nal and transverse structure functions {olIgt)- It is re- 
lated to the other structure functions by. 



R{x,Q^) 



F2 iM'x 



1 



2xFi'' Q2 ' - 2xFi 
where Fl is called the longitudinal structure function, 

UPx^ 



(2) 



Tl{x,Q')^T2 1 + 



Other useful relations are: 



2xFi. 



2xT^(x, Q2) = Q2^ 1 + ^M^^^IQ'' 



l + i?(a;,Q2) ' 



(3) 



(4) 



l + i?(a;,Q2)- 



The three structure functions T^ix.Q'^^, Fi{x,Q'^) and 
xT^(x,Q'^) depend on x and . 

Integrating over x, the differential dependence on u 
can be written in the simplified form 



da"''' _ I B V C v"^ \ 
dv y^^lfy^7L2El)' 



(5) 



The coefficients A, B, and C depend on integrals over 
structure functions, where 



Jo 

B = J [T2{x) T xT3{x)yx, 

GlM 



(6) 



C = B- ^^-^ I T2{x) R dx, 
'o 



TT 



and 



R 



1 + ^ Mx 



l + R{x,Q^) V 



In the limit vjEi, — )• 0, the A term dominates and the 
B and G terms are very small. The MINOS collaboration 
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used the number of "low-z^" events (with v < Vcut) in the 
detector to determine the relative flux of neutrinos and 
antineutrinos as a function of . 

In the MINOS analysis, the relative flux is determined 
using events with v < 1 GeV for energies in the range 
3 < < 9 GeV, and for in the range 5 < i?^ < 9 GeV. 
Events with v < 2 GeV are used for i'^ and events in 
the range 9 < -Ejy < 18 GeV, and events with v < 5 GeV 
are used for > 18 GeV. 

MINOS divides the number of "low-z^" events with y < 
Ucut = Vcut/Ei, by correction term fc to account for the 
energy dependence from the B and C terms. Here 



1.05 



/c(i?.) 



C 

A2EI 



dy. (7) 



As seen in equation |6j the negative contribution of J-2{x) 
in B partially cancels the positive contribution of xJ-^{x) 
for t'^'s. For D^^s both contributions are negative. There 
are additional small corrections that are applied to equa- 
tion [6] to correct for differences in IF2 between neutrinos 
and antineutrinos. 

In practice, a neutrino interaction generator model [H] 
is used to compute fc from 



jcyi^u) — —, — ^ r- 



(8) 



and 



The corrections factors fc used by MINOS for 
are shown in Fig. [T] 

The measured "low-z/" sample is corrected for detector 
smearing and acceptance by multiplying the number of 
observed "low-i^" events in the data in each energy bin by 



R^'^^iE,) 



which is obtained from a Monte Carlo detector simula- 
tion. JV^BAf jyREC g^j,g ^]^g number of generated and 
reconstructed events below Vcut in each reconstructed en- 
ergy bin, respectively. In the first pass the initial input 
flux from a beam model is used. It is then replaced by the 
extracted "low-i^" flux and the procedure is reiterated to 
account for the effect of the flux model on the acceptance 
corrections. (The change in the extracted flux is found to 
be negligible) . The "low-z^" sample is further corrected for 
radiative effects using Ref. [3]. The absolute level of the 
flux is set by normalizing the cross section in data to a 
nominal world average charged current cross section at 
some high energy. As mentioned earlier, in MINOS the 
normalization is set to the average of previous (Jtotai/E 
measurements for neutrino energies between 30 and 50 
GeV. 

There are three criteria for the effectiveness of the 
"low-j/" method. 

1. The number of "low-i/" events that are used in the 
determination of the flux should not be a large fraction 
of the total number of neutrino events in each energy 
bin. 
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Fig. 1. The "low-;/" correction factors fc used by MINOS for 
neutrinos (shown in the top panel) and antineutrinos (shown 
in the bottom panel). 



2. The systematic uncertainty in the energy dependent 
correction factor fc should be small. 

3. The number of "low-z^" events that are used in the 
determination of the flux should be sufficiently large 
to have flux sample with small statistical errors. 

The first two criteria require a Vcut which is as low as 
possible. The third requires a Vcut which is as large as 
possible. 

The MINOS collaboration uses the criteria that the 
fractional contribution of events with v < Vcut to the total 
charged current cross section should be less than 60%. 
MINOS uses events with v <1 GeV for determination of 
the flux at their lowest f^, energies. 

The fraction of events with v < 1 GeV is less than 
60% for interactions with i?^ > 3 GeV and for 
interactions with Ep > 5 GeV. Therefore, to determine 
the flux for E^ < 3 GeV and i?p < 5 GeV we need to use 
a Vcut which is smaller than 1 GeV. 

We investigate Vcut = 0.25 GeV to be used for i^i^.p > 
0.7 GeV, and 

i^cut — 0.5 GeV to be used for Ei, p > 1.4 
GeV. These samples can be cross calibrated against the 
i^cut — 1 GeV sample in the range E^ > 3 GeV for neu- 
trinos and Ep > 5 GeV for antineutrinos. Similarly, they 
can be calibrated against the Vcut = 2 GeV and v^ut = 5 
GeV samples in the range E^^i> > 9 GeV and E^^p > 18 
GeV, respectively. 
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Fig. 2. The kinematic region in the (in GeV^), u = Ehad (in 
GeV) plane for v^^ (/J^) energies less than 4 GeV. The shaded 
area \s v < 0.25 GeV (color online). 
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Fig. 3. The kinematic region in the Q2 (in GeV^), 
(in GeV) plane for i/^ {Vfj) energies less than 1.5 GeV. The 
shaded area \s v < 0.25 GeV (color online). 



kinematic region in (in GeV^ and v = Ehad (in GeV) 
for E^<4: GeV. 

Fig. [s] shows the kinematic region in and for 
El, < 1.5 GeV. In this paper we focus on i/ <0.25 GeV 
region (shaded area in Fig. [2] and Fig. |3]) for the lowest 
neutrino energies. In addition, we investigate the <0.50 
GeV region as an additional check. 

For E^= 3 GeV, about 1/3 of the total charged cur- 
rent cross section originates from QE scattering, 1/3 from 
resonance production and 1/3 from inelastic scattering. 

As seen in Fig. |3] the < 0.25 GeV sample is domi- 
nated almost entirely by QE events with < 0.45 GeV^. 

The v < 0.5 GeV sample includes both QE events with 

< 0.95 GeV^ and also Z\(1232) resonance events with 

< 0.3 GeV^. Both samples include a very small fraction 
of events originating from coherent pion production (as 
discussed in Appendix II). 

In the very "low-z^" region it is more convenient to 
write the expression for the charged current differential 
cross sections as follows [IOllTT| : 



dQ^dv 



Sc 



(1 



'2E^ 



m,, 



AE^ 



2ME AE ME 



f 4M2£:2 

,2 



-ScosW5 



m. 



ME 



(9) 



where S^os = ^ cos^ 0c = 80 x 10"^^° crnVGeV^. In 
the scattering process, there are additional small contri- 
butions from strangeness and charm non-conserving pro- 
cesses. In the discussion below we do not show these terms 
explicitly, but charm and strangeness changing contribu- 
tions are assumed to be included in the analysis. (The 
strangeness changing valence quark contributions are pro- 
portional to ^ sin^ 9c) ■ 

Each of the structure functions has a vector and axial 
component (except for Wy, which originates from axial- 
vector interference). The vector part of Wi and are 
well known since they are related to the vector part of W2 
and Wi by the following expressions [lUj: 



3 The "\o\M-v" method at low energies 

In the few GeV region, there are several types of neutrino 
interaction processes as defined by the final state invariant 
mass W . These include quasielastic (QE) reactions [W < 
1.07 GeV), production of the Zi(1232) resonance (1.1 < 
< 1.4 GeV), coherent pion production, production of 
higher mass resonances (1.4 < W < 2.0 GcV) and the 
inelastic continuum {W > 2.0 GeV ). Fig. [2] shows the 



y^vector -y^vector 



y^vector y^vector 



M'v' 



Mv 



.M' 



Q 



2 ■ 



At 'low-z/' and very high energy the charged current 
cross section is only a function of 1^2 ■ If we integrate the 
cross section from Vmin ~ up to 1^= Vcut (where Vcut is 
small), we can write the expression for the cross section in 
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terms of only, and energy dependent corrections ratios 
to the component: 



(fa 



dQ^dv 



(10) 



(74 + 0-5, 



Here, 



^ (JW2 (00), where 

(TW2 = Sec 

(JW2 (00) = Sec 



Wo dv. 



W2 dv. 



and the small corrections to the QE cross section are 

Q2 + r 



O2 — Se 



<Jl—Se 



= Se 



0-4 — Sc, 



CT5 = •S'c, 



V 

E 



W2 dv. 



(11) 
(12) 

(13) 



Wi dv. 



2ME AE ME 

Wa dv 



W3 dv. 



ME 



W5 di/. 



The above can be written in terms of fractional correc- 
tions: 



o-i^cut{E) = <TW2{'^) [fc] , 

fc = [fw2 + /2 + /l +f3 + h + h] 
<^W2 



(14) 



71^2 = 



aw2 (00) 
(TW2 (00) 



1, 



The energy dependent corrections fw2, /i, /2, fs, /s, 
and fi and /s can be calculated within a specific models. 
The theoretical uncertainty in fc determines the system- 
atic uncertainty in the relative flux which can extracted 
from the "low-i^" events. 



1- fw2 — 



1 is well known and does not con- 



1JW2 (00) 

tribute to the uncertainty in fc- 

2. The energy dependent correction /2 is explicit and 
therefore does not contribute to the uncertainty in fc- 

3. The contributions of fi and /s are small since they 
are proportional to the square of the muon mass, and 
therefore have a negligible contribution to the uncer- 
tainty in fc. (Note that the vector parts of /4 and 
/5 are known very well since they can be expressed in 
terms of the vector parts of Wi and W2)- 

4. The only non-negligible uncertainty originates from 
the modeling of the contributions of /i and /a (pri- 
marily from /a). 
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Fig. 4. Top panel: Neutrino partial charged cross sections per 
nucleon for "low-v" events determined from the GENIE Monte 
Carlo [12]. Also shown are the measurements of MINOS on 
iron (per nucleon corrected for the excess number of neutrons). 
Bottom panel: The fraction of "low-i^" neutrino events in the 
GENIE 12 Monte Carlo as compared with the measurements 
of MINOS' (color online). 



The technique does not depend on the modeling of W2 
because the aw2 cross section is the same at all energies. 
All energy dependent corrections are expressed in terms of 
ratios to aw2 ■ In quark parton language, the uncertainty 
in /i is related to the uncertainty in the longitudinal struc- 
ture function at low and the uncertainty in /a is related 
to the uncertainty in level of antiquarks in the nucleon at 
low Q^. For QE scattering and resonance production the 
structure functions are expressed in terms of form factors. 



3.1 Partial charged current cross sections 

The top panel of Fig. [4]shows the partial neutrino charged 
current cross section per nucleon for "low-i^" events (for v 
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Same as Fig. [4] for the case of antineutrinos. (color 



cuts of 0.25, 0.5, 1, 2 and 5 GeV) as a function of energy as 
determined by the GENIE[:12 Monte Carlo for a carbon 
target. The top panel of Fig. [5] shows the corresponding 
partial charged current cross sections for antineutrinos. 

Also shown are the measurements of the partial charged 
current cross sections on iron from the MINOS collabora- 
tion (for 1/ cuts of 1, 2 and 5 GeV). The MINOS cross 
sections for iron have been corrected for the excess num- 
ber of neutrons in iron. Note that the nuclear corrections 
to the structure functions in iron nucleus are larger than 
in carbon. Therefore, the partial cross sections on carbon 
and on iron may not be the same. 

At high energies (as shown in Fig. 4]and[5]) the partial 
cross sections for a fixed Vcut are independent of energy 
and are approximately equal for neutrinos and antineu- 
trino. The fact that these partial charged current cross 
section are relatively independent of energy is the basis 
for the "low-z/" method. 



The bottom panels of figures |4] and [5] show the fraction 
of "low-!/" events predicted by the GENIE Monte Carlo 
as compared with the measurements in MINOS. In order 
to use the technique at low energies the fractions must 
be smaller than 0.6. Therefore, at the lowest energies we 
must use v cuts of 0.25 and 0.50 GeV. 

MINOS is a sampling target calorimeter which has 
poor resolution at low hadron energy. Therefore, "low-i/" 
samples with h' < 0.25 GeV and v < 0.5 GeV cannot be 
defined reliably. On the other hand, since the MINERvA 
detector is a fully active target calorimeter, "low-:^" sam- 
ples with u < 0.25 GeV and v < 0.5 GeV can be used. 



3.2 Absolute normalization 

Since the neutrino energy range for MINERvA is limited 
to lower energies, we propose that the MINERvA charged 
current cross section measurements be normalized to the 
cross section in the energy range between 10 to 20 GeV 
(e.g. at a mean energy of 15.1 GeV). The absolute level of 
the charged current cross section at this energy range has 
been measured by both the MINOS and NOMAD collab- 
orations. 

The MINOS total cross section measurement for an 
isoscalar iron target at a neutrino energy of 15.1 GEV is 

fjMiNos 1^ = 0.708 ± 0.020 X lO^^^cmVCeV 

per nucleon in iron. Here the total error of 0.02 is the 
combined statistical, systematic and normalization errors 
of 0.008 ± 0.012 ± 0.015, respectively. 

The NOMAD cross section measurement for an isoscalar 
carbon target at a neutrino energy of 15.1 GEV is 



NOMAD 



IE = 0.698 ± 0.025 x lO^^'^cmVGeV 

per nucleon in carbon. 

The MINOS total cross section measurement for an 
isoscalar iron target at an antineutrino energy of 15.1 GEV 
is 

^MiNOS 1^ = 0.304 ± 0.012 X lO^^^cmVCeV 

per nucleon in iron. Here, the total error of 0.012 is the 
combined statistical, systematic and normalization errors 
of 0.007 ± 0.007 ± 0.006, respectively 

Alternatively, it may be possible for MINERvA to nor- 
malize to the partial cross sections measured by MINOS 
for i/ < 1 GeV and v <2 GeV at 15.1 GeV. These partial 
cross sections (which were used by MINOS to determine 
their relative flux) are relatively constant between 10 and 
20 GeV. However, the MINOS partial cross sections are 
measured on iron. The MINERvA target is solid scintil- 
lator (i.e. carbon), and the partial cross sections for iron 
and carbon can be different. For a neutrino energy of 15.1 
GeV MINOS measured the following isoscalar partial cross 
sections on iron (per nucleon): 



MINOS 



(15.1) = 1.729 ± 0.049 x 10~-^**cm^(z^ < 2 GeV) 
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Fig. 6. The v < 0.25 GeV partial charged current cross sec- 
tions (per nucleon) as a function of energy from the GENIE 
Monte Carlo (for carbon target). Shown are the QE contri- 
bution, the contribution from pion production process (e.g. 
A, inelastic and coherent pion) and the total. The u < 0.25 
GeV cross sections for i/^ are shown on the top panel, and the 
f < 0.25 GeV cross sections for are shown on the bottom 
panel, (color online). 

^^"^^■^(IS.l) = 0.968 ± 0.027 x 10"^**cm2(i/ < 1 GeV). 

For an antineutrino energy of 15.1 GeV MINOS has 
measured the following isoscalar partial cross sections: 



^MINOS 



(15.1) = 1.585 ± 0.063 x lQ~^^cw?[v < 2 GeV) 
M/7voS(-j^5j-) = 0.939 ± 0.039 x 10-^^cm^{i^ < 1 GeV). 



and the ly < 0.25 GeV cross sections for are shown 
on the bottom panel. The QE contribution is shown in 
red, the contribution from pion production process (e.g. 
A, inelastic and coherent pion) is shown in blue and the 
total is shown in black. Most of the events are QE and the 
contribution from pion production processes is negligible. 

As mentioned earlier, the technique does not rely on 
the modeling of >V2, or the modeling of nuclear effects 
(e.g. Fermi motion smearing) on W2. This is because the 
cross section (including nuclear effects) is the same 
at all neutrino energies. 

The uncertainty in the flux extracted from the event 
sample with i/ < 0.25 GeV is determined by how well we 
can model the relative contributions of Wi and for the 
case of QE scattering on bound nucleons, or equivalently 
the relative contributions of /i and to fc ■ Here /i and 
/3 are proportional to the ratios ^ and Since the 

ratios ^ and ^ for QE scattering on free nucleons are 
very well known, the uncertainty in fc originates primar- 



ily from modeling the nuclear corrections to ^ and 
for nucleons bound in a nuclear target. 
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Fig. 7. The transverse enhancement ratio[T3] (7?.t) as a func- 
tion of . Here, TZt is ratio of the integrated transverse re- 
sponse function for QE electron scattering on nucleons bound 
in carbon divided by the integrated response function for in- 
dependent nucleons. The black points are extracted from Carl- 
son et aiUni, and the blue bands are extracted from a fit|19| 
to QE data from the JUPITER[T8] experiment (Jlab exper- 
iment E04-001). The curve is a fit to the data of the form 
Ht = 1 + AQ^e-'^ The dashed lines are the upper and 
lower error bands (color online). 



4 Using "\o\M-u" events with u < 0.25 GeV 

As seen in Fig. |3]the v < 0.25 GeV region is dominated 
by QE events. This is illustrated in Fig. [6] which shows 
the relative contributions of QE and non-QE processes 
to v < 0.25 GeV cross section as a function of energy 
(as determined from the GENIE Monte Carlo). The ly < 
0.25 GeV cross sections for are shown on the top panel. 



4.1 Quasielastic t'^,^'^ scattering 

The relationship between the structure functions and form 
factors for v^,v^ QE scattering[13j on free nucleons is 
given bv P3T5] : 
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Fig. 8. Comparison of predictions for the total QE 

cross section sections from the nominal TE model, the "In- 
dependent Nucleon (MA=1.014)" model, the "Larger Ma 
(Ma=1.3) model", and the "QE+np-nh RPA" MEC model 
of Martini et al.|^ The data points are the measurements of 
MiniBooNE[2T] (gray stars) and NOMAD[22] (purple circles) 
(color online). 
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Fig. 9. Contribution of the various components {o'W2, ""2, cri, 
(T3, (74, (75 ) to the total QE cross section (as predicted by the 
TE model). Top panel: Neutrinos. Bottom panel: Antineutrinos 
(color online). 
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where 



an(i 



r[gXim? 

l+T 



Here, G^(Q2), GHQ^), GIj{Q^) and G'^V(Q') are the 
electric and magnetic nucleon form factors, which are mea- 
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sured in electron scattering experiments. Note that: 

(Jrp 



^'^'"^ oc r|gr,(g2)|2; a-^"'^' cx (1 + t)\Ta{Q^)? 



^axial Q 



Therefore, for QE i/f^, scattering only Q'^j contributes 
to the transverse virtual boson absorption cross section. 



4.2 Transverse enhancement QE scattering from nuclei 

Studies of QE electron scattering on nuclear targets[l6] 
indicate that only the longitudinal part of the QE cross 
section can be described in terms of a universal response 
function of independent nucleons bound in a nuclear po- 
tential (and free nucleon form factors). In contrast, a sig- 
nificant additional enhancement with respect to the model 
is observed in the transverse part of the QE cross section. 

The enhancement in the transverse QE cross section 
has been attributed [TB] to meson exchange currents (MEC) 
in a nucleus. Within models of meson exchange currents [IS] 
the enhancement is primarily in the transverse part of the 
QE cross section, while the enhancement in the longitu- 
dinal QE cross section is small (in agreement with the 
electron scattering experimental data). 

The conserved vector current hypothesis (CVC) im- 
plies that the corresponding vector structure function for 
the QE cross section in v^, scattering can be expressed 
in terms of the structure functions measured in electron 
scattering on nuclear targets. Therefore, there should also 
be a transverse enhancement in neutrino scattering. In 
models of meson exchange currents the enhancement in 
the axial part of v^i^^v^ QE cross section on nuclear tar- 
gets is also expected small. 

The transverse enhancement observed in electron scat- 
tering is a function of both and v. However, a simple 
way to account for the integrated transverse enhancement 
from nuclear effects is to assume that G\j{Q'^) and G'^j{Q'^) 
are enhanced in a nuclear targets by factor ^/Rtl- 

Bodek, Budd and Christy [11] have used electron scat- 
tering data[16 y i"8 l fT9] to parametrize Rtl as follows: 



with A = 6.0 and B = 0.34 GeV^. The electron scat- 
tering data indicates that the transverse enhancement is 
maximal near (5^=0.3 GeV^ and is small for greater 
than 1.5 GeV^. The upper error band is given by A = 6.7 
and B = 0.35 GeV^, and the lower error band is given by 
A = 5.3 and B = 0.33 GeV^. 

In modeling j/^ QE scattering on nuclear targets we 
use BBBA20Q725 parameterization ,17J of the free nucleon 
electromagnetic form factors G^((5^), C^iQ"^), G\j{Q'^) 
and GljiQ'^) (with = 0.71 GEV^), and a dipole axial 
form factor with Ma — 1.014 GeV. We apply the trans- 
verse enhancement correction to G\j{Q^) and G\j{Q^). 
We also apply Pauli blocking corrections to the differential 
QE cross section as parametrized by Paschos and Yujill^. 
We refer to this model as the Transverse Enhancement 



(TE) model. This is the nominal model that is used in 
this paper. 

We also compare calculations based on the nominal 
TE model to two other models. The first model is the in- 
dependent nucleon model with Pauli blocking with Ma = 
1.014 GeV, without transverse enhancement. We refer to 
this model as the "Independent Nucleon (MA=1.014)" 
model. This model, which is used by the NOMAD 
collaboration, is very close to the model which is cur- 
rently implemented in the GENIE Monte Carlo (the GE- 
NIE defauh value is Ma = 0.99 GeV). The second model 
is the independent nucleon model with Pauli blocking. 
Ma = 1.3 GeV, without transverse enhancement. This 
model is used by the MiniBooNE Collaboration[2T]. We 
refer to this model as the "Larger Ma (Ma=1.3) model". 
We use the difference between the three models as a con- 
servative systematic error on the flux extracted from the 
v samples. 

Fig. [8] shows a comparison of predictions of various 
model predictions for the u^, total QE cross section sec- 
tions to experimental data on nuclear targets. Shown are 
"Independent Nucleon (MA=1.014)" model, the "Larger 
Ma (Ma=1.3) model", and the TE model (with upper 
and lower error bands). Also shown are the predictions of 
the "QE-hnp-nh RPA" MEC model of Martini et al.[2n] 
The data points are the QE cross section measurements 
of MiniBooNE [21 (gray stars) and NOMAD _22^ (purple 
circles). Note that there is an overall « 10% systematic 
error in the experimental QE cross sections because of 
uncertainties in the determination of the neutrino and an- 
tineutrino fluxes in each of the two experiments. 

In this paper we use the error band in the transverse 
enhancement parameters as a lower limit on systematic 
error in the modeling. We use the "Independent Nucleon 
(MA=1.014)" and the "Larger Ma (Ma=1.3) model" as 
conservative upper limits on the errors in the modeling. 

Fig. [9| shows the contribution of the various compo- 
Inents {crw2^ '^2, cri, crs, cr4, CT5) to the total QE cross 
section (as defined by Eq. 14) as a function of incident 
energy. These contributions are calculated using the TE 
model. The top panel shows the contribution of the var- 
ious components for the neutrino QE cross section, and 
the bottom panel shows the contribution of the various 
components for the antineutrino QE cross section. 



4.3 Neutrino and antineutrino total cross sections 

The MINOS collaboration uses the criteria that the frac- 
tion of "low-i^" events that are used for the determination 
of the relative neutrino flux in an energy bin should be less 
than 60% of the total number of charged current events. In 
order to test for this fraction, we need to use a parameter- 
ization to estimate the energy dependence of the neutrino 
and antineutrino charged current total cross sections. 

10 and 11 show the i'^ and total charged cur- 



Fig. 



rent cross sections measured on isoscalar nuclear targets 
by the MINOS |3j (iron), NOMAD [22J (carbon), and Serpukov]: 
(Serp96, aluminum) experiments. The total cross sections 
per nucleon (divided by neutrino energy) are shown in 
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Fig. 10. 

BNL82[21 

on isoscalar nuclear targets for v in units of 10~'^'*cm'^/GeV. 
The orange line shows the predictions of the unmodified GE- 
NIE Monte Carlo. The QE cross section in the GENIE MC is 
shown as the blue line. The QE contribution calculated with 
the TE model is shown as a green line. The dashed blue line 
shows the prediction of the modified GENIE MC (using the 
TE model QE cross section instead). The thick brown line is a 
parameterization described in the text (color online). 



units of 10~'^*cm^/GeV (with statistical, systematic and 
normalization errors combined in quadrature). The ratio 
of the and z/^ total charged current cross sections is 
shown in the bottom panel of Fig. [TT] The cross sections 
reported by the MINOS collaboration were measured us- 
ing a neutrino flux extracted from "low-i^" samples with 
V less than 1, 3, and 5 GeV. 

T0| are low energy cross sections 
(BNL82). Since the BNL82 cross 
sections were measured on a deuterium target we apply 
a correction to account for nuclear effects. The BNL82 
points shown in the figure were increased by the differ- 
ence of the predictions of the TE model for the QE cross 
section (which is expected to describe the cross section 
on a heavy nuclear target) and the "Independent Nucleon 
(MA=1.014)" model (which is expected to describe the 
QE cross sections on deuterium). 

The orange line shows the predictions of the GENIE 
Monte Carlo. The QE cross sections in the GENIE MC 
are computed using the independent nucleon model with 
Ma = 0.99 GeV. The QE contribution to the cross section 
from GENIE is shown as a blue line. The QE contribu- 
tion calculated with the TE model is shown as a green 
line. The curve labeled GENIE with QE-TE (shown as 
a dotted blue line) represents the GENIE cross section 
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measured by at BNLj^] 
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Fig. 11. Same as Fig. 
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but (a) for the antineutrino charged 
current cross section, (b) for the ratio of antineutrino and neu- 
trino total cross sections (color online). 



increased by the difference of the predictions of the TE 
model for the QE cross section (which is expected to de- 
scribe the cross section on a heavy nuclear target) and 
the "Independent Nucleon (MA=0.99)" model (which is 
currently implemented in GENIE). 

In our investigation of the "low-z^" technique, we use a 
parameterization to estimate the total i^^, charged cur- 
rent cross sections. The parameterization, which is shown 
as the thick red line in Fig. flO) is given by 



^ = [A + B e-^"'^^+D 
E,, 



where for we use A^= 0.675, 0.12, Cl,.= 9 GeV, 
=0.4, C2^= 3 GeV2, C3^= 0.22 GeV, and K = 1.0. 
For Df, we use ^^,=0.329, Bp= -0.06 and Clp=I3 GeV 
Di,^ 0.09, C2p= 30 GeV2 , C3^^ 0.8 GeV, and K = 0.8. 
Here, ^ is total charged current cross section per nucleon 

in units of 10~^*cm2/GeV. 
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Fig. 12. Contribution of the various components {pwi, ""2, o\, 
az, (74, 0-5) to the v < 0.25 GeV partial charged current cross 
section. This sample is dominated by QE i'fj,N — > M'f' events. 
Top panel: Neutrinos. Bottom panel: Antineutrinos (color on- 
line) . 



Tlie above form is constrained to yield the average 
world cross section measurements in the 30 to 50 GeV re- 
gion of of 0.675 10"38pj^2/QgY^ 0.329 lO-^ScmVGeV 
for and i^^, respectively. 

We only use this parameterization to estimate the frac- 
tional contribution of "low-z/" events to the total cross 
section to determine the region where it is less than 60%. 
When improved total cross section measurements become 
available (e.g. from MINERvA), this parameterization can 
be updated to include the new data. 



4.4 Results with v < 0.25 GeV 



Fig. 12 shows the contribution of the various components 
{<Jw2^ ^2, ci, (T3, (74, (75) to the V < 0.25 GeV partial cross 
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Fig. 13. The v < 0.25 GeV sample for z/^ scattering on carbon. 
Top panel: The total correction factor fc (black line), the con- 
tribution of the kinematic correction to W2 (/2) (yellow line), 
the contributions from Wi (/i) (pink line), the contribution 
from W3 (fs) (blue line), and the very small contributions of 
VV4 ifi), and W5 (/s). Bottom panel: The fractional contribu- 
tion of 1/ < 0.25 GeV events to the total 1/^ charged current 
cross section, (color online). 



section. This sample is dominated by QE z^^V — )• fi^ P 
events. The partial cross section as a function of energy for 
neutrinos is shown in the top panel and the partial cross 
section for antineutrinos is shown in the bottom panel. 
The partial cross section (per nucleon) is calculated on a 
carbon target using the TE model. 

The uncertainty in the relative values of the v < 0.25 
GeV partial cross section as a function of energy deter- 
mines the uncertainty in the determination of the relative 
fluxes. Here fc{E) is the ratio of the partial cross section 
to the value of the partial cross section at E — 00. 

Fig. 13 ^a) (top) shows the correction factor fc for the 



V < 0.25 GeV sample for neutrinos as a function energy. 
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Fig. 14. Same as Fig. 13 for the case of antineutrinos. 
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Fig. 15. The error band in the correction factor /c for v < 
0.25 GeV. Top panel: Neutrinos. Bottom panel: Antineutrinos 
(color online). 



The error bands in fc (originating from the uncertainty 
in the transverse enhancement) are shown as the dotted 
lines, and represent the lower limit on errors. Also shown 
on the figure is the negative contribution from the kine- 
matic correction /2 (which is well known), and the contri- 
butions of /i, /s, /4 and f^. Here the contribution of 
and /s is negligible. For the case of neutrino scattering, 
the positive contributions of /i and fs partially cancel the 
negative contribution of /2. Fig. 13 b) (bottom) shows the 
fractional contribution of the < 0.25 GeV sample to the 
total neutrino charged current cross section. This fraction 
is less than 60% for z/^ energies above 0.70 GeV. 

Fig. [14] is the same as Fig. [13] for the case of antineu- 
trinos. For the case of antineutrino scattering /a changes 
sign, and both /2 and /a are negative. The fractional con- 
tribution of the 1/ < 0.25 GeV sample to the total antineu- 
trino charged current cross section is less than 60% for 
energies above 1.0 GeV. 



4.4.1 Uncertainty in the fc correction factors 

It has been traditional to use the value and error in the 
effective Ma extracted from neutrino scattering data as an 
estimate of various uncertainties. Typically, the difference 
between resuhs with Ma = 1.014 GeV and Ma = 1.3 
GeV are used an upper limit on the error. 

We find that the values of the fc correction factor are 
insensitive to Ma- This is because at small Q^, both ratios 
/i, and /a are insensitive to AIa- Specifically, both 
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Fig. 16. Comparisons of our calculated values of the nor- 
malized /c:>.<o.25(15.1)(£;) (=/c(15.1) for u < 0.25 GeV) to 
values from the GENIE MC. The values calculated with the 
nominal TE model for QE scattering {Ma = 1.014 GeV) are 
shown in black. The values calculated assuming no transverse 
enhancement and AIa ~ 1.014 GeV are shown in red. The 
GENIE prediction (which has no transverse enhancement and 
uses AIa = 0.99 GeV) is close to the red curve as expected, 
(color online). 



are insensitive to because the change in at small 
is small. Since fc is insensitive to large variations in 
AIa one may naively surmise that the error in fc is small. 

However, we find that the difference between the val- 
ues fc calculated with and without transverse enhance- 
ment is larger than the error estimate extracted from the 



is sensitive to 



uncertainty in AIa- This is because 

OXiiQ^): which depends on the magnitude of the trans- 
verse enhancement at small Q^. 

Fig. 15 shows the errors in fc from the uncertainty 



tainties in the TE parameters is also very small (less than 
0.005). 

We obtain a more conservative estimate of the system- 
atic error in fc originating from uncertainties in the mod- 
eling the QE cross section by taking the difference between 
fc calculated with and without transverse enhancement. 
At the lowest energy of 0.7 GeV, this difference is -0.05 
for i/^. Since at 0.7 GeV fc ~ 1.3 this corresponds to 
a maximum error in the determination of the flux of 
3.8%. 

For the difference between fc calculated with and 
without transverse enhancement at an energy of 1.0 GeV 
is -1-0.03. Since at 1.0 GeV fc ~ 0.6 this corresponds a 
maximum error in the determination of the flux of 5%. 



4.5 Comparison to GENIE and fc. 



(15.1 GeV) 



We have used a sample of events generated by the GENIE 
Monte Carlo. Our studies are done at the generated level 
and therefore do not depend on the detector parameters 
or energy resolutions of any specific experiment. 

We extract the energy dependence of the v < 0.25 
GeV cross section from the GENIE MC sample using the 
following expression: 



'i^<0.25 



where the superscript AIC refers to events generated by 
the GENIE Monte Carlo. 

Here, N^'^^{E) is the number of events generated by 
the Monte Carlo with neutrino energy E, and N^'^ 25{E) 
is the subset of these events with ly < 0.25 GeV. 

As mentioned earlier, we propose that the neutrino 
cross sections at low energy be measured relative to the 
neutrino cross section at 15.1 GeV. For any cross section 
model we can define the normalized quantity fc v<o 25(15-1 
GeV) as: 



/c:.<o.25(15.1)(F) = 
which is equivalent to 

/C:.<0.25(15.1)(F) 



T.<o.25(S= 15.1 GeV) 

fC:u<0.2b{E) 

fc:.<o.25{E = 15.1 GeV) 



in the TE parameters. The error originating from uncer- 



We compare the values of /c:i^<o. 25(15. l)(i?) predicted by 
the GENIE MC to our calculations. 

For completeness, we give the values of /c:i/<o. 25(15.1) 
that can be used to convert between fc:v<Q.2b{E) and 

fC:u<0.25{E). 

For the TE QE model we find /c:i.<o.25(15.1)=-1.018 
(for v) and 0.966 (for v). For QE models without TE we 
find similar values of /c:i.<o.25(15.1)=1.016 (Ma=1.014) 
and /c:..<o.25(15.1)=1.014 (Ma=1.03) for v. For v we find 
/c:^<o.25(15.1)=0.969 for models without TE. 

Comparisons of our calculated values of the normal- 
ized /c:i.<o.25(15.1)(£') to values from the GENIE MC 
are shown in Fig.|16| The top panel shows the comparison 
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for neutrinos and the bottom panel shows the compari- 
son for antineutrinos. Our calculation for the TE model 
is shown in black. Our calculation assuming no transverse 
enhancement and Ma = 1.014 GeV is shown in red. As 
mentioned earlier, the values for Ma = 1.014 GeV (red 
line) and Ma = 1.3 GeV (blue line) are very close to each 
other. The GENIE prediction (which has no transverse 
enhancement and uses Ma = 0.99 GeV) is close to the 
red curve as expected. The GENIE predictions include a 
contribution from coherent pion production. As shown in 
Appendix II, for the ly < 0.25 GeV sample, the contribu- 
tion from coherent pion production is less than 0.1% for 
neutrinos and less than 0.6% for antineutrinos. 



4.6 Conclusions of the studies with v <0.25 GeV 

In conclusion, we find that the method works very well for 
v < 0.25 GeV. If one takes the average of all the models, 
a conservative upper limit of the model uncertainty in the 
relative flux extracted from the f < 0.25 GeV sample 
is 1.9% for energies above 0.7 GeV and 2.5% for i/^ 
energies above 1.0 GeV. The GENIE Monte Carlo is in 
reasonable agreement with the models and therefore can 
be used to obtain a first order neutrino flux. 

A study of the distributions of QE events in MIN- 
ERvA can be used to constrain the dependence of the 
QE differential cross sections and thus reduce the model 
dependence in the determination of the relative flux to a 
negligible level. A GENIE Monte Carlo which is tuned to 
agree with the new data can be used to extend the tech- 
nique to lower energies. 

5 Using "low-z/" events with u < 0.5 GeV 

The 1/ < 0.5 GeV and Df^ samples have close to twice 
the number of events as the v < 0.25 GeV samples. These 
samples for scattering are also dominated by QE events, 
but include a significant fraction (about 1/3) of events 
in which a single pion is produced in the final state. As 
seen in Fig. |3] the v < 0.5 GeV samples are composed of 
QE events with < 0.9 GeV^ and Z\(1232) events with 
Q2 < 0.3 GeV2. 

Fig. [TT] shows the i' < 0.5 GeV partial charged current 
cross sections as a function of energy. The partial cross sec- 
tions extracted from the GENIE Monte Carlo are shown 
as black points with MC statistical errors. The v < 0.5 
GeV partial cross section for i/^ scattering is shown on 
the top panel, and the ly < 0.5 GeV partial cross section 
for scattering is shown on the bottom panel. The QE 
contribution to the < 0.5 GeV partial cross section is 
shown in red, and the contribution from pion production 
processes (Z\, inelastic and coherent pion production) is 
shown in blue. 

As seen in Fig. |17[ the pion production contribution to 
the < 0.5 GeV partial cross section is relatively constant 
with energy, while the QE contribution has some energy 
dependence. Therefore, the energy dependence of the sum 
of the two contributions to the v < 0.5 GeV partial cross 
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Fig. 17. The v < 0.5 GeV partial charged current cross section 
as a function of energy from the GENIE Monte Carlo. The 
QE contribution is shown in red, the contribution from pion 
production process {A, inelastic and coherent pion production) 
is shown in blue, and the total u < 0.5 GeV GeV partial cross 
section is shown in black. The f < 0.5 GeV GeV partial cross 
section for i/^i is shown in the top panel, and the i/ < 0.5 GeV 
partial cross section for i^^ is shown in the bottom panel, (color 
online). 



section requires modeling of the relative magnitude of QE 
and pion production processes (specifically at low Q^). 
As shown in Fig. 



18 and [19] the consistency among 



the experimental measurements of pion production cross 
sections in the region of the Z\(1232) resonance is about 
20% (depending on the neutrino energy and the nuclear 
target). We use this variation to get an estimate of the 
model uncertainty in the determination of the neutrino 
flux from the v < 0.5 GeV samples. This uncertainty can 
be greatly reduced when more precise measurements of the 
QE and pion production cross sections become available 
(e.g. from MINERvA). 
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Fig. 18. u^P II Zi++ (top panel) and u^N ii+ A 
(bottom panel) cross sections (for W < lA GeV) measured 
on free nucleons (H and D), compared to predictions from the 
GENIE MC (black points with errors), (color online). 



6 Pion production with W < 1.4 GeV 

In this section we describe the uncertainties in the mod- 
eling of pion production cross sections for W < lA GeV. 
The antineutrino structure functions are related to the 
neutrino structure functions by the following relationship. 



(15) 



•'i •'i 

6.1 i/^P ^ H-A++ and i?^N /i+Z\- (FIT-A) 



We define the cross section for i/^P — > /i~Z\++ as the 
integrated cross section for W < lA GeV for the following 
single final state: 
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Fig. 19. Same as Fig, 
online). 



Paschos2011(M^=1.05,c"=1.2) 
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for the case of nuclear targets, (color 



We define the cross section for D^N — ^ fj,^A~ as the in- 
tegrated cross section for W < 1.4 GeV for the following 
single final state: 



-Ntt- 



Therefore, our definition includes the sum of the con- 
tributions of the resonant cross section and the non-resonant 
continuum. 

The structure functions (form factors) for the reac- 
tions Vfj^P — ^ fi~A~^'^ and i/^TV — >• defined above 
are the same (except that for antineutrinos the structure 
function W3 changes sign). It has been experimentally 
determined j32, that f^P cross section for W < 1.4 GeV is 
dominated by the resonant production process. Simi- 
larly, the < 1.4 GeV cross section for P^^N is dominated 
by the resonant A~ production process. 

As discussed in the Appendix, we parametrize the Z\+"'" 
and A~ production cross sections in terms of form factors 
as given by Paschos and Lalakulich|10|. with the form fac- 
tors of Paschos and Schalla[10] . In order to obtain predic- 
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Fig. 20. v^,N fi A+ (top panel) and v^,P -> /i+Zi" (bot- 
tom panel) cross sections [W < 1.4 GeV) measured on free 
nucleons (H or D). The predictions from the GENIE MC are 
shown as black points with errors, (color online). 



tions for the W < lA GeV region, we divide all theoreti- 
cal A production cross sections by a factor of 1.2 (because 
20% of the resonant cross section is above W = lA GeV). 
We vary two of the parameters in the model, specifically 
and to obtain a band that span the experimental 
data. We extract from the measured distributions 
and use to set the overall normalization. 



The top panel in Fig. [18] shows a summary of cross 
section measurements for v^P — > /i^Z\++ on free nu- 
cleons (hydrogen or deuterium targets.) Shown are bub- 
ble chamber measurements at low energy from Argonne 
(ANL73[1S], ANL79[11], ANL82[27]) and measurement at 
low energy from Brookhaven (BNL86[28|). Also shown are 
measurements at higher energies from the Fcrmilab bub- 
ble chamber (FNAL78[29J, FNAL81 30 ) and high energy 
data from CERN (BEBC80 31J, BEBC80[32^). The bot- 
tom panel in Fig. [l8| shows the BEBC90|32] cross section 
measurements for P^A^ — >■ fj,~^A~ on free nucleons (deu- 



5 

E 

a 45 



3.5 



x10-'» 



2.5 



1.5 



0.5 







-FIT-B(M =1.62,c*= 

A 5 


1.27] 


GENIE Carbon 


















































































































*i — \ — — 




















/ 1 




























+ 
























































































i 









6x10' 1 2 3 4 5 6 7 8 910 

vp A Cross Section W<1 .4 on Carbon 



20 30 40 50 
E(GeV) 




2 3 4 5 6 7 8 910 



20 30 40 50 
E(GeV) 



Fig. 21. Same as Fig. 20 for nuclear targets, (color online). 



terium target). The predictions from the GENIE MC on 
free nucleons (shown as black points with MC statistical 
errors) are near the upper bound of our three parameter- 
izations. 

The black curve labeled Paschos-2011 (Mf =1.05, 
— 1.2) uses the original values of AI^ and from the 
paper |10j by Paschos and Lalakulich. These values were 
obtained from fits to cross sections and distributions 
measured at low energies at Brookhaven and Argonne. 
The red curve labeled FIT-Al (Mf =1.93, = 0.62) is 
derived from a fit to the cross sections and distribution 
of the higher energy BEBC90^ data for Df,N -> n+A-. 
The blue curve labeled FIT-A2 (Mf =1.75, = 0.49) is 
derived from a fit to the cross sections and distribution 
of the higher energy BEBC90[32 data for i/^P fi- A++ . 

The top panel in Fig. [19] shows a summary of cross 
section measurements for v^P — > fj,~ A'^'^ data on nuclear 
targets. Shown are the measurements of Gargamelle78 [55] 
(Propane), SKAT88[51| (Freon), and SKAT89(55] (Freon). 
The bottom panel shows measurements of h'i_iN — > /i"'"Z\~ 
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cross sections on nuclear targets from Gargamelle78[36] 
(Propane) and SKAT89[3F (Freon). 

Aside from Pauli suppression and final state interac- 
tion, the structure functions (form factors) for the pro- 
cesses in Fig. [Ts] and 19 are the same. The black (Paschos- 



V A Cross Section per Nucleon W<1 .4 on Carbon 



2011)^ed (FIT-Al) and blue (FIT-A2) curves shown in 
Fig. \m and \m use the free nucleon form factors (but in- 
clude the Pauli suppression for the case of nuclear tar- 
gets). The calculations do not include the effect of final 
state interaction for the nuclear targets. The three curves 
(Paschos-2011, FIT-Al and FIT-A2) conservatively span 
all the available A'^'^ and A~ production cross sections 
on hydrogen, deuterium and nuclear targets, as shown in 



Fig. [18] and 19 The cross sections for the production of 
Z\"'"+ and on nuclear targets predicted by GENIE are 
near the upper bound of our three parameterizations. Ad- 
ditional details are given in the Appendix. 



6.2 i,^N 



'Z\+ and D^P 



^I+A° (FIT-B) 



We define the cross section for ly^j^N — /i^Zi+ as the sum 
of the integrated cross sections for W < lA GcV for the 
following two final states: 



fl-NTT+ 

' frPn" 



We define the cross section for i/^P 



/i+Z\*^ as the sum 
of the integrated cross sections for W < lA GeV for the 
following two final states: 



Therefore, our definition includes the sum of the con- 
tributions of the resonant cross section and non-resonant 
continuum. 

The structure functions (form factors) for the reac- 
tions i^fj_N — >• fi~A~^ and D^j^P — >■ n~^A° defined above 
are the same (except that for antineutrinos the structure 
function W3 changes sign). Because of Clebsch-Gordan 
coefficients [TIT the resonant cross section for Z\+ produc- 
tion in i^f^N collisions is a third of the resonant cross sec- 
tion for Z\++ production in i/^P collisions. Similarly, the 
resonant cross section for A^ production in D^P collisions 
is a third of the cross section for resonant production of 
A~ in Df^N collisions. 

However, unlike the case for f^P {A'^^) and i^^iV 
{A~), where the cross sections are dominated by the res- 
onant process, there is a significant contribution from the 
non-resonant continuum to the W < lA GeV cross section 
in I'p^N and D^^P collisions. 



The top panel of Fig. 



20 



shows the ly^N 



cross sections {W < 1.4 GeV) measured on free nucleons 
(deuterium). Shown are measurements from ANL79[26]. 
ANL82[23, and BEBC90i32j. The predictions from the 
GENIE MC are shown as black points with MC statisti- 
cal errors. In order to describe the data (which has a large 
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Fig. 22. The total cross sections on carbon (per nucleon) 
predicted by GENIE for W < lA GeV (black points with MC 
statistical errors) for f^C — >■ {fi~ A'^'^ or A^) are shown on 
the top panel, and for i/^C — >■ fi'^{A° or A~) are shown on the 
bottom panel, (color online). 



non-resonant contribution) we changed the parameters in 
the Paschos and Lalakulich^lO resonance model to fit the 
observed distribution and total W < lA GeV cross 
sections. The green curve labeled FIT-B (Mf =1.62, 
~ 1.27) is derived from a fit to the W < 1.4 GeV cross 
sections and distribution of the BEBC90[32 data for 
i>p,N — >■ /i~Z\+. This curve provides a parameterization 
which describe the experimental data for the production 
of A'^ (with neutrinos) and Z\° (for antineutrinos) on free 
nucleons. The GENIE MC cross sections for the produc- 
tion of Z\+ on free nucleons are lower than the fit. 



The structure functions (form factors) for the reactions 
i^^N fi-A+ and j>^P fi+A° {W < 1.4 GeV) are are 
same. The bottom panel of Fig. [20| shows a comparison of 
the predictions of FIT-B (Mf =1.62, = 1.27) (green 
curve) for the i?^P ^i^ A^ cross sections on free nucleons 
compared to the predictions from the GENIE MC which 
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Fig. 23. The v < 0.5 GeV sample for v^. This sample includes 
both QE VfiN — >■ fi~ P events («66%) and A production events 
(« 33%). Top panel: The total corrections factor fc (with error 
bands) and the contributions of the kinematic correction to 
W2 (/2), and the contributions from Wi (/i), W3 (/a), W4 
(/4), and W5 (/s). Bottom panel: The fractional contribution 
of < 0.5 GeV events to the total cross section, (color online). 
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Fig. 24. Same as Fig. 23 for the case of antineutrinos. (color 
online). 



6.3 Comparisons of < 1.4 GeV cross sections on 
carbon 



are shown as black points with MC statistical errors. The 
GENIE MC cross sections for the production of on free 
nucleons are lower than the fit. 



Fig. 



20 



shows the prediction of FIT-B (Mf =1.62, 
— 1.27)"(g'reen curve) for the v^N — >• ^~zi+ (top panel) 
and Pfj,P — ?► /^+zi° (bottom panel) W < 1.4 GeV cross 
sections on nuclear targets compared to predictions from 
the GENIE MC (black points with MC statistical errors). 
The cross sections on nuclear targets are expected to be 
somewhat lower than the cross sections on free nucleons 



(which are shown in Fig. 20 1 . Here, FIT-B includes the ef- 
fect of Pauli suppression (but not final state interaction). 
The GENIE MC cross sections for the production of Z\+ 
and A'^ on nuclear targets are lower than the fit. Addi- 
tional details are given in the Appendix. 



A more relevant comparison is to determine how well the 
GENIE Monte Carlo describes the sum of the proton and 
neutron cross sections on carbon, since it is the total num- 
ber of < 0.5 GeV events on carbon that are used in the 
determination of the neutrino flux. 

Fig. [22] shows the predictions from the GENIE MC 
for total A production cross section for W < lA GeV 
on carbon (per nucleon). The neutrino cross sections for 
I'^C — ^ /i^(Z\++ or Z\+) are shown in the top panel, and 
the antineutrino cross sections i^^C — )■ (^"*"Z\° or A^) are 
shown in the bottom panel. The cross sections which are 
predicted by GENIE arc compared to our three parame- 
terizations. (Paschos-2011, FIT-Al and FIT-A2 for A++ 
and A-, and FIT-B for A+ and A°). The GENIE cross 
section predictions for the total A production cross sec- 
tions on carbon (which use the Rein and Seghal model [37] 
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Fig. 25. Comparisons of our calculated values of the normal- 
ized /c:>.<o.5(15.1)(£) (=/c(15.1) for V < 0.5 GeV) to val- 
ues from the GENIE MC. Our nominal model (shown as the 
solid black line) uses the TE model for QE scattering and the 
Paschos 2011 model for A production. Neutrinos are shown 
on the top panel and antineutrinos are shown on the bottom 
panel (color online). 
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Fig. 26. The error band in the normalized correction factor 
/c:;.<o.5(15.1)(£;) (=/c(15.1) for u < 0.5 GeV). Our nominal 
model is QE with transverse enhancement and the Paschos 
2011 model for A production. Shown are the differences be- 
tween our nominal model and other model assumptions for 
neutrinos (top panel) and for antineutrinos (bottom panel), 
(color online). 



for resonance production) fall near the lower bound of our 
three parameterizations of the experimental data. 

As described below, the uncertainties in the measure- 
ments of the A production cross sections do not place a 
serious limitation on the flux extractions using the low v 
method. 



6.4 Determination of neutrino and antineutrino flux 
using J/ < 0.5 GeV samples on carbon 

The ly < 0.5 GeV sample includes both QE v^N — > fi^ P 
events (w 66%) and A production events (ss 33%). 

The top panel of Fig. [23] shows the total correction 
factor fc{E) for the v < 0.5 GeV sample (defined as 
fc:v<Q.b{E)) for neutrino running. Also shown are the 



various contributions to fc-.vKO.siE) including the kine- 
matic correction to (/2), and the contributions from 
Wi (/i), W3 (/s), W4 (/4), and W5 (/5). The bottom 
panel shows the fractional contribution oi ly < 0.5 GeV 
events to the charged current neutrino total cross section. 
Using our nominal model (TE model for QE scattering 
and the Paschos 2011 model for A production) we find 
that the fraction of < 0.5 GeV events is less than 60% 
for energies above 1.2 GeV. 

The top panel of Fig. [24] shows the total correction fac- 
tor fc:v<o.5{E) for antineutrino running. Also shown are 
the various contributions to fc:v<o.b including the kine- 
matic correction to 'W2 (/2), and the contributions from 
Wi (/i), W3 (/a), W4 (A), and W5 (/s). The bottom panel 
shows the fractional contribution of < 0.5 GeV events to 
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the charged current antineutrino total cross section. Using 
our nominal model (TE model for QE scattering and the 
Paschos 2011 model for A production) we find that the 
fraction of j/ < 0.5 GeV events is less than 60% for 
energies above 2 GeV. 

As for the ly < 0.25 sample, we propose that the neu- 
trino and antineutrino cross sections at low energy be mea- 
sured relative to the cross sections at 15.1 GeV. Therefore, 
we define normalized quantity fc:v<o.b 

{l^.l){E) for the 

V < 0.5 sample as: 

/c:.<o.5(15.1)(£;) = a,<QAE)/a,<:o.5{E = 15.1 GeV) 
which is equivalent to 

/c:.<o.5(15.1)(^) = fc{E)/fciE = 15.1 GeV) 

Here, /c:z.<o.5(E=15.1 GeV)=1.0113 (for i^) and 0_.9507 
(for P). These values can be used to convert between fcM<o.5 
and fc:v<o.b{E). 

Fig. [25] show^ comparisons of our calculated values of 
the normalized fc:u<o.b{^^-^){E) (shown as the solid black 
line) to values extracted from the GENIE MC. The GE- 
NIE predictions include a contribution from coherent pion 
production. As shown in Appendix II, for the u < 0.50 
GeV sample, the contribution from coherent pion produc- 
tion is less than 0.7% for neutrinos and less than 3% for 
antineutrino 

Our values are calculated from our nominal model which 
uses the TE model for QE scattering and the Paschos 2011 
model for A production. Neutrinos are shown on the top 
panel and antineutrinos are shown on the bottom panel 
(color online). 

Figure [26] shows the error band in the correction factor 
/c:y<o.5(15.1)(i?) for neutrinos (top panel) and antineu- 
trinos (bottom panel). The error band is defined as the 
differences between our nominal model and other model 
assumptions. For neutrinos with energies greater than 1.2 
GeV, the error in /c(15.1) is less than 0.03, which corre- 
sponds to a 2.6% upper limit on the model uncertainty in 
the neutrino flux extracted from the < 0.5 GeV sample. 
For antineutrinos with energies greater than 2 GeV the 
error in /c(15.1) is less than 0.01 (which corresponds to 
a 1.4% upper limit on the model uncertainty in the an- 
tineutrino flux extracted from the < 0.5 GeV sample). 

In order to go to lower neutrino and antineutrino ener- 
gies we need to use the i' < 0.25 GeV sample. The model 
uncertainty in the relative flux extracted from the < 0.25 
GeV sample is 1.9% for energies above 0.7 GeV and 
2.5% for Pfj, energies above 1.0 GeV. With improved de- 
termination of QE and A production cross sections (e.g. 
in MINER vA), the model uncertainties can be further re- 
duced, and the method may be extended to lower energies. 



6.5 Resolution, acceptance and radiative corrections 

The i/ < 0.25 GeV events are primarily QE events with 
Q2 < 2M X 0.25 « 0.45 GeV^. We can select either all 
events with ly < 0.25 GeV or only QE events with < 0.5 
GeV2 
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Fig. 27. da/dQ^ cross sections (for W < 1.4 GeV) measured 
on deuterium at high energies by AUasia et al. fBEBC90|32)V 
The cross sections for for j/^P — > A^'^ are shown on the top 
panel and the cross sections for i>f_iN — > n"^ A~ are shown on 
the bottom panel, (color online). 



The ratio of the number of reconstructed events with 
ly < 0.25 GeV (or < 0.5 GeV^) in data and MC as 
a function of energy is proportional to the ratio of the 
true flux to the simulated flux in the MC. This ratio pro- 
vides a measure of the relative neutrino flux as a func- 
tion of energy. A complete Monte Carlo should include the 
small contributions from coherent pion production (shown 
in Appendix II), strange particle production such as QE 
production of hyperons[38j, and radiative corrections [391 
HP] . The effects of experimental resolution and acceptance 
should also be simulated. 

At present the GENIE Monte Carlo includes coherent 
pion production, but does not include the QE production 
of hyperons, nor radiative effects. 

If the GENIE Monte Carlo is used, then one may wish 
to weight the rate of QE events (as a function of Q^) by the 
ratio of events expected in the TE model to the number 
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Fig. 28. Top panel: da/dQ^ for —> /i^/i^^ cross sections 
(for W < 1.4 GeV) measured on hydrogen at high energies by 
Allen et. al. (BEBC80[3T]). Bottom panel: da/dQ^ for Uf^N 
fj,~ cross sections (for W < 1.4 GeV) measured by BEBC90 
on free nucleons on deuterium, (color online). 



of events predicted by the model which is implemented 
in GENIE (i.e. the "Independent Nucleon" model with 
Ma — 0.99 GeV). In addition, QE production of hyperons 
and radiative effects need to be added. 



7 Conclusions 

We find that the model uncertainties in using the "low-i^" 
event samples with ly < 0.25 and ly < 0.5 GeV are well 
under control (less than 3%). Therefore, the "low-i/" tech- 
nique can be used at low energies (0.7 GeV for neutrinos 
and 1 GeV for antineutrinos). Once data from MINErVA 
on QE scattering and resonance production becomes avail- 
able, the model uncertainties can be made even smaller, 
and the technique may be extended to even lower energies. 

Since the model uncertainties are under control, the 
dominant systematic error originates from how well the 



detector response is understood. Specifically, the mis- re- 
construction of high v events as "low-z/" events must be 
modeled reliably. This is because at high energies (as shown 
in Fig. [2| mis-reconstruction of the hadron energy of high 
v events can increase the number of "low-i/" events, while 
at low energies there are fewer high u events that can be 
mis-reconstructed at low v. 

The dominant uncertainty in the method comes from 
the calibration and resolution smearing in the measure- 
ment of the hadronic energy. This was the dominant er- 
ror when this method was used in MINOS because of the 
poor resolution of the MINOS target calorimeter at low 
hadronic energy. 

As mentioned in the introduction, the standard method 
for the determination of the neutrino flux requires the 
modeling of pion production as well as the complicated 
magnetic focusing elements. The determination of the flux 
for the Fermilab NUMI beam with the standard method is 
limited at present by the uncertainties in pion production 
cross sections. The resulting error in the flux is about 5% 
at low energies (1-2 GeV) and 10%-15% at the higher ener- 
gies (10-20 GeV). Therefore, having the "low i/" method 
which yields the relative neutrino flux as a function of 
energy very useful. In principle, the uncertainties in the 
standard method can be improved with better measure- 
ments of pion and kaon production cross sections. Plans 
for such future measurements at the GERN Laboratory 
are currently under discussion. 

A second method, which requires the measurement of 
the muon rate downstream of the decay pipe, can not de- 
termine the energy dependence of the flux. It mostly con- 
strains the overall level of the flux. At present, the uncer- 
tainties in the overall calibration of the muon chambers 
yields an uncertainty in the flux of about 10%. . 

A third method uses inverse muon decay + e — > 
li~ -^Vf, events in the detector. The threshold for this reac- 
tion is about 12 GeV. Therefore, this method can only be 
used at higher energies. Inverse muon decay was used by 
NOMAD to constrain their neutrino flux at high energies. 
In addition to being statistically limited, the final state 
energy of inverse muon decay events is not fully measured 
since there is a neutrino in the final state. This places a 
limitation on the determination of the energy dependence 
of the neutrino fluxes. This method cannot be used for the 
determination of the flux for antineutrinos. 

A fourth method uses the neutral current reaction v^^ -\- 
e — ^ z/^ -|- e. In addition to being statistically limited, the 



final state energy in 



i^u + e events is not fully 



measured since there is a neutrino in the final state. This 
places a limitation on the determination of the energy 
dependence of the fluxes. In this method only the sum 
of the fluxes for neutrinos and antineutrinos is measured 
because calorimetric detectors such as MINERvA cannot 
determine the charge of final state electrons. 
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Fig. 29. The three v^P/vi_^N da/dQ^ cross sections models 
(for W < 1.4 GeV) with Pauli suppression for nuclear targets 
at an energy of 40.5 GeV. The cross sections for for i/^P — > 
H~ A'^^ are shown on the top panel and the cross sections for 
i>fiN — )■ jj.'^A~ are shown on the bottom panel, (color online). 



8 Appendix I: A production cross sections 
8.1 A production form factors 

For the vector contribution we use the formulae for the 
structure functions Wi, W2, W3, W4, and W5 on free nu- 
cleons from Lalakulich and Paschos|lO]. We neglect the 
effect of Fermi motion. The form factors that we use are 
taken from Paschos and Schalla^Oj. Specifically, the vec- 
tor form factors are 
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Fig. 30. Our da/dQ^ cross sections model (for W < 1.4 GeV) 
with Pauli suppression for nuclear targets at an energy of 40.5 
GeV. The cross sections for for v^N — >■ ^'^ A~ are shown on the 
top panel and the cross sections for S-^P — >■ /i^zi" are shown on 
the bottom panel. The green curve labeled FIT-B (Ma=1.62, 
Ci = 1.27) represents a fit to the BEBC90 Vf,N free nucleon 
data (color online). 

with My = 0.84 GeV, which have been extracted from 
electroproduction data. 

For the vector-axial interference Wz{Q^,v) Paschos 
and Schalla use the form factor C^{Q'^) where 



1 



(1 



Here, we define = C^(0) 

Paschos and Schalla use low energy 7r+p — > where 
the non-resonant background is smallest. With M^=1.05 
GeV they extract value of C^(0) = 1.08 from the data. 
Since this value is close to 1.20 predicted by the Goldberger- 
Treiman relation, they chose to use Ci^=1.2. 
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Paschos and Schalla mention that several recent arti- 
cles also calculate (0) by fitting experimental data [¥T1 
[42 .43,44. 45 ■„46..47, with values varying from 0.87 up to 
1.20. Models with a resonant background [41,42 prefer 
the power value, while the other articles P3i,44i,45.i4^IT7] 
prefer values closer to 1.20. The reasons for the differences 
is the treatment of the non-resonant background, the form 
of the axial form factor that is used, and the exact kine- 
matics at small . 

For Z\++ and A~ we define the Pachos-2011 param- 
eterization using above form factors with C^=1.2 (ex- 
tracted through PCAC), M^=1.05 GeV, and the vector 
form factors described above. As mentioned earlier, FIT- 
Al and FIT-A2 use the same form but with different val- 
ues of and Ma- 

For Z\+ and Z\° production our Fit-B uses the same 
form factors multiplied by a factor of l/\/(3) (as expected 
from Clebsch-Gordan coefficients [TU]). However, in order 
to account for the large non-resonance background, we use 
different values and Ma- 



8.2 Various parameterizations 

The low energy and high energy data for neutrino and 
antineutrino production of the Z\(1232) resonance are not 
entirely consistent. Therefore, we use range of parameter- 
ization to span the systematic error in our modeling of A 
production cross sections. 

The form factors for v^P /i" 
IX^ A~ should be the same. The da/dQ'^ differential cross 
sections {W < 1.4 GeV) for u^P — >• A~^^ measured at 
high energies are shown in the top panel of Fig. [27] (Al- 
lasia et. al., BEBC90,32j data on deuterium) and also on 
the top panel of Fig. [2| (Allen et. al. BEBC80[HT] data on 
hydrogen) The bottom panel of Fig. 27 shows the da/dQ"^ 



CC Coherent Pion Fraction on Carbon for v (v<0.25) 



'A++ and i>^,N 



cross sections at high energies {W < 1.4 GeV) for v^^N — ^ 
/i+zl^ measured by AUasia et. al. (BEBC90) data on deu- 
The black curve labeled Paschos-2011(Mf =1.05, 
1.2) is from fits to lower energy f^P — > /i^Z\++ 
data (BNL and Argonne). The red curve labeled FIT-Al 
(Ma=1.93, = 0.62) is a fit to the BEBC90 P^TV ^ 
^+A- data . The blue curve labeled FIT-A2 (AfA=1.75, 
= 0.49) is a fit to the BEBC90 v^,P ^ A++ data. 
The variation among the three curves is taken as a sys- 
tematic error. 

The bottom panel of Fig. 28 shows values of da/dQ^ 
differential cross sections for Vi^i_N — > A^ (for W < 
1.4 GeV) measured by BEBC90 on free nucleons on deu- 
terium. This reaction has different form factors then v^P — > 
H~A++- The green curve labeled FIT-B (Ma=1.62, 
= 1.27) represents a fit to the BEBC90 i/f^N data. 

We use the above models with the addition of Pauli 
suppression in order to model the differential cross sections 
on nucl ear targets. 

shows our three v^P/v^N da/dQ^ cross sec- 



Fig. 
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tions moHels (for W < l.A GeV) with Pauli suppression 
for nuclear targets at an energy of 40.5 GeV. The cross 
sections for for v^P — > /i^Z\++ are shown on the top 
panel and the cross sections for D^N — > fj,~^ A~ are shown 
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Fig. 31. The fraction of events from coherent pion production 
in the u < 0.25 GeV event sample (calculated with GENIE) 
as a function of neutrino energy: Neutrinos (top panel) and 
antineutrinos (bottom panel). 



on the bottom panel. These two reactions should be de- 
scribed by the same form factors. The black curve labeled 
Paschos-2011(Mf =1.05, = 1.2) is from fits to lower 
energy f^P free nucleon data (BNL and Argonne). The 
red curve labeled FIT-Al (Ma=1.93, = 0.62) is from 
a fit to the BEBC90 D^j^N free nucleon data. The blue curve 



labeled FIT-A2 {Ma^I-75. 



0.49) is from a fit to the 



BEBC90 I'fj.P free nucleon data. The variation among the 
three curves is taken as a systematic error. 

Fig. pffl shows our v^N/v^P dajdQ^ cross sections 
model (for W < 1.4 GeV) with Pauli suppression for nu- 
clear targets at an energy of 40.5 GeV. The cross sections 
for for i/^iV — > /i+zi^ are shown on the top panel and the 
cross sections for v^P — >■ /z+zi" are shown on the bottom 
panel. The green curve labeled FIT-B (A/4=1.62, 



1.27) is extracted fit to the BEBC90 
data. 



,N free nucleon 
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Fig. 32. Same as 31 but for the ly < 0.50 GeV sample. 



9 Appendix II: Coherent Pion Production 

Fig. [3l] shows the fraction of events from coherent pion 
production in the ly < 0.25 GeV event sample (calculated 
with GENIE) as a function of neutrino energy. Neutrinos 
are shown on the top panel and antineutrinos are shown 
on the bottom panel. For the ly < 0.25 GeV sample, the 
contribution from coherent pion production is less than 
0.1% for neutrinos and less than 0.6% for antineutrinos 

Fig. [32] shows the fraction of events from coherent pion 
production in the ly < 0.50 GeV event sample (calculated 
with GENIE) as a function of neutrino energy. Neutrinos 
are shown on the top panel and antineutrinos are shown 
on the bottom panel. For the ly < 0.50 GeV sample, the 
contribution from coherent pion production is less than 
0.7% for neutrinos and less than 3% for antineutrinos 
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